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This letter describes how to perform searches over the complete kinematically-allowed parameter 
space for new pair-produced color octet particles that each subsequently decay into two jets plus 
missing energy at the Tevatron. This letter shows that current searches can miss otherwise discov- 
erable spectra of particles due to CMSSM-motivated cuts. Optimizing the Ht and JEt cuts expands 
the sensitivity of these searches. 



I. INTRODUCTION 

In many theories beyond the Standard Model, there 
is a new color octet particle that decays into jets plus a 
stable neutral singlet. This occurs, for example, in su- 
persymmetry [1] and Universal Extra Dimensions [2], as 
well as Randall-Sundrum 3 and Little Higgs models [4 . 
As a result, jets plus missing transverse energy (fix) is 
a promising experimental signature for new phenomena 

nana® ei . 

At present, the jets + fix searches at the Fermi- 
lab Tevatron are based upon the minimal supersymmet- 
ric standard model (MSSM) and look for production of 
gluinos (g) and squarks (q), the supersymmetric partners 
of gluons and quarks, respectively [7J [5]. Both gluinos 
and squarks can decay to jets and a bino (B), the su- 
persymmetric partner of the photon. The bino is stable, 
protected by a discrete R-parity, and is manifest as miss- 
ing energy in the detector. Different jet topologies are 
expected, depending on the relative masses of the gluinos 
and squarks. 

There are many parameters in the MSSM and set- 
ting mass bounds in a multidimensional parameter space 
is difficult. This has lead to a great simplifying ansatz 
known as the CMSSM (or mSUGRA) parameterization 
of supersymmetry breaking |10j . This ansatz sets all 
the gaugino masses equal at the grand unified scale and 
runs them down to the weak scale, resulting in an ap- 
proximately constant ratio between the gluino and bino 
masses (m^ : rrig = 6 : 1). Thus, the mass ratio be- 
tween the gluino and bino is never scanned when search- 
ing through CMSSM parameter space. Since the bino is 
the LSP in most of the CMSSM parameter space, the 
restriction to unified gaugino masses means that there is 
a large region of kinematically-accessible gluinos where 
there are no known limits. 

The CMSSM parametrization is not representative 
of all supersymmetric models. Other methods of super- 
symmetry breaking lead to different low-energy particle 
spectra. In anomaly mediation [11 , the wino can be 



the LSP; for instance, rrig : ra^ ~ 9:1. Mirage 
mediation [IS], in contrast, has nearly degenerate gaug- 
inos. A more comprehensive search strategy should be 
sensitive to all values of rrig and rrig. Currently, the 
tightest model-independent bound on gluinos is 51 GeV 
and comes from thrust data at ALEPH and OPAL [13.. 

In this paper, we describe how bounds can be placed 
on all kinematically-allowed gluino and bino 1 masses. 
We will treat the gluino as the first new colored particle 
and will assume that it only decays to the stable bino: 
g — > qiq* — > q\qiB. The spin of the new color octet 
and singlet is not known a priori; the only selection rule 
we impose is that the two have the same statistics. In 
practice, the spin dependence is a rescaling of the entire 
production cross section. For our analysis, we will as- 
sume that the octet has spin 1/2, and will show how the 
results vary with cross section rescaling. 

We show how a set of optimized cuts for fix and 
Ht = Sjets can discover particles where the current 
Tevatron searches would not. In order to show this, we 
model our searches on D0 's searches for monojets [5], 
squarks and gluinos [7J. In keeping the searches closely 
tied to existing searches, we hope that our projected sen- 
sitivity is close to what is achievable and not swamped 
by unforeseen backgrounds. 



II. EVENT GENERATION 
A. Signal 

The number of jets expected as a result of gluino 
production at the Tevatron depends on the relative mass 
difference between the gluino and bino, rrig — m^. When 
the mass splitting is much larger than the bino mass, 



Throughout this note, we will call the color octet a "gluino" 
and the neutral singlet the "bino," though nothing more than 
the color and charge is denoted by these names. 
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FIG. 1: Boosted gluinos that are degenerate with the bino 
do not enhance the missing transverse energy when there is 
no hard initial- or final-state radiation. (A) illustrates the 
cancellation of the bino's pT- (B) shows how initial- or final- 
state radiation leads to a large amount of fir even if the 
gluino is degenerate with the bino. 

the search is not limited by phase space and four or 
more well-separated jets are produced, as well as large 
missing transverse energy. The situation is very differ- 
ent for light gluinos (nig < 200 GeV) that are nearly 
degenerate with the bino. Such light gluinos can be co- 
piously produced at the Tevatron, with cross sections 
O(10 2 pb), as compared to 0(1O~ 2 pb) for their heav- 
ier counterparts (nig > 400 GeV). Despite their large 
production cross sections, these events are challenging 
to detect because the jets from the decay are soft, with 
modest amounts of missing transverse energy. Even if 
the gluinos are strongly boosted, the sum of the bino 
momenta will approximately cancel when reconstruct- 
ing the missing transverse energy (Fig. To discover 
a gluino degenerate with a bino, it is necessary to look 
at events where the gluino pair is boosted by the emis- 
sion of hard QCD jets (Fig. |TJ3) . Therefore, initial-state 
radiation (ISR) and final-state radiation (FSR) must be 
properly accounted for. 

The correct inclusion of ISR/FSR with parton show- 
ering requires generating gluino events with matrix ele- 
ments. We used MadGraph/MadEvent [2] to compute 
processes of the form 

PP^ 99 + Nj, (!) 

where N = 0,1,2 is the multiplicity of QCD jets. The 
decay of the gluino into a bino plus a quark and an anti- 
quark, as well as parton showering and hadronization of 
the final-state partons, was done in PYTHIA 6.4 [15]. 

To ensure that no double counting of events occurs 
between the matrix-clement multi-parton events and the 
parton showers, a version of the MLM matching proce- 
dure was used [IB] . In this procedure, the matrix el- 
ement multi-parton events and the parton showers are 
constrained to occupy different kinematical regions, sep- 
arated using the k± jet measure: 

d 2 (i,j) = Ai^-min^,^) 
d 2 (i,beam) = p 2 ^, (2) 



where ARfj = 2(coshA7y - cosA^) [T7]. Matrix- 
element events are generated with some minimum cut- 
off d(i,j) = Q^in- After showering, the partons are 
clustered into jets using the kr jet algorithm with a 
Qmfn > Qmin- The event is then discarded unless all re- 
sulting jets are matched to partons in the matrix-element 
event, <i(parton, jet) < Qj£f n . For events from the high- 
est multiplicity sample, extra jets softer than the soft- 
est matrix-element parton are allowed. This procedure 
avoids double-counting jets, and results in continuous 
and smooth differential distributions for all jet observ- 
ables. 

The matching parameters (<9^ and Q^fn) should 
be chosen resonably far below the factorization scale of 
the process. For gluino production, the parameters were: 

QME = 20 GeV and Q^f n = 30 GeV. (3) 

The simulations were done using the CTEQ6L1 
PDF [W\ and with the renormalization and factoriza- 
tion scales set to the gluino mass. The cross sections 
were rescaled to the next-to-leading-order (NLO) cross 
sections obtained using Prospino 2.0 [T§] . 

Finally, we used PGS 20J for detector simulation, 
with a cone jet algorithm with AR = 0.5. As a check 
on this procedure, we compared our results to the signal 
point given in '?! and found that they agreed to within 
10%. 



B. Backgrounds 

The three dominant Standard Model backgrounds 
that contribute to the jets plus missing energy searches 
are: W ± /Z° + jets, ti, and QCD. There are several 
smaller sources of missing energy that include single top 
and di-boson production, but these make up a very small 
fraction of the background and are not included in this 
study. 

The W ± /Z° + nj and ti backgrounds were gen- 
erated using MadGraph/MadEvent and then showered 
and hadronized using PYTHIA. PGS was used to recon- 
struct the jets. MLM matching was applied up to three 
jets for the W /Z° background, with the parameters 

QME = 1Q Q eV and gP^ = 15 Ge y Xhe tQp back _ 

ground was matched up to two jets with = 14 GeV 
and Qjnf n = 20 GeV. Events containing isolated leptons 
with pt > 10 GeV were vetoed to reduce background 
contributions from leptonically decaying W bosons. To 
reject cases of JUt from jet energy mismeasurement, a 
lower bound of 90° and 50° was placed on the azimuthal 
angle between $t and the first and second hardest jets, 
respectively. An acoplanarity cut of < 165° was applied 
to the two hardest jets. Because the DQ analysis did 
not veto hadronically decaying tau leptons, all taus were 
treated as jets in this study. 
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Simulation of the missing energy background from 
QCD is beyond the scope of PYTHIA and PGS, and was 
therefore not done in this work. However, to avoid the 
regions where jet and calorimeter mismeasurements be- 
come the dominant background, a lower limit of J£t > 
100 GeV was imposed. Additionally, in the dijet analy- 
sis, the azimuthal angle between the Pt and any jet with 
Pt > 15 GeV and \rj\ < 2.5 was bounded from below by 
40°. This cut was not placed on the threejet or multijet 
samples because of the large jet multiplicities in these 
cases. 

For each of the W /Z° + nj and tt backgrounds, 
500K events were generated. The results reproduce the 
shape and scale of the fix and Ht distributions pub- 
lished by the D0 collaboration in |7J for lfb -1 . For 
the dijet case, where the most statistics are available, 
the correspondence with the D0 result is ±20%. With 
the threejet and multijet cuts, the result for the tt back- 
ground is similar, while the W ± /Z° + nj backgrounds 
reproduce the D0 result to within 30 — 40% for the 
threejet and multijet cases. The increased uncertainty 
may result from insufficient statistics to fully populate 
the tails of the $t and Ht distributions. The PGS prob- 
ability of losing a lepton may also contribute to the rel- 
ative uncertainties for the W ± + nj background. Heavy 
flavor jet contributions were found to contribute 2% to 
the W ± /Z° backgrounds, which is well below the uncer- 
tainties that arise from not having NLO calculations for 
these processes and from using PGS. 



III. PROJECTED REACH OF SEARCHES 

A gluino search should have broad acceptances over 
a wide range of kincmatical parameter space; it should 
be sensitive to cases where the gluino and bino are nearly 
degenerate, as well as cases where the gluino is far heav- 
ier than the bino. As already discussed, the number of 
jets and Pt depend strongly on the mass differerence be- 
tween the gluino and bino. Because the signal changes 
dramatically as the masses of the gluino and bino are 
varied, it is necessary to design searches that are gen- 
eral, but not closely tied to the kinematics. We divided 
events into four mutually exclusive searches for Pt plus 
lj, 2j, 3j and 4 + j, respectively. For convenience, we 
keep the nj + Pt classification fixed for all gluino and 
bino masses (see Table |l| . These selection criteria were 
modeled after those used in D0 's existing search [7]. 2 
These exclusive searches can be statistically combined to 
provide stronger constraints. 
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TABLE I: Summary of the selection criteria for the four non- 
overlapping searches. The two hardest jets are required to be 
central (|jy| < 0.8). All other jets must have \n\ < 2.5. 



Two cuts are placed on each search: H T ™ and 
$r mm . In the D0 analysis, the H T and p T cuts are 
constant for each search. The signal (as a function of 
the gluino and bino masses) and Standard Model back- 
ground are very sensitive to these cuts. To maximize the 
discovery potential, these two cuts should be optimized 
for all gluino and bino masses. For a given gluino and 
bino mass, the significance (S/s/ S + B) is maximized 
over Hf in and ^ T min in each nj + p T search. Due to 
the uncertainty in the background calculations, the S/B 
was not allowed to drop beneath the conservative limit of 
S/B > 1. More aggressive bounds on S/B may also be 
considered; D0 , for instance, claims a systematic uncer- 
tainty of O(30%) in their background measurements [7]- 
The resulting 95% sensitivity plot using the optimized 
Ht and Pt cuts is shown in Fig. [2] The corresponding 
inset illustrates the effect of varying the production cross 
section. 

For light and degenerate gluinos, the lj + Pt and 
2j + Pt searches both have good sensitivity. In an inter- 
mediate region, the 2j + p Tl 3j + Pt and 4+j + p T all 
cover with some success, but there appears to be a cov- 
erage gap where no search does particularly well. If one 
does not impose a S/B requirement, a lot of the gap can 
be covered, but background calculations are probably 
not sufficiently precise to probe small S/B. For massive, 
non-degenerate gluinos, the 3j + Pt and 4+j + p T both 
give good sensitivity, with the 4 + j + Pt giving slightly 
larger statistical significance. 

In the exclusion plot, the Pt and Ht cuts were op- 
timized for each point in gluino-bino parameter space. 
However, for gluino masses 200 GeV % nig < 350 GeV, 
where the monojet search gives no contribution, we 
found that the exclusion region does not markedly 
change if the following set of generic cuts are placed: 



2 It should be noted, however, that the D)0 searches are inclu- 
sive because each is designed to look for separate gluino/squark 
production modes (i.e., pp — > qq, qg, gg). 



(Ht,Pt) > (15O,lOO) 2i+ T , 

(150, 100) 3]+fiT , (200, 100) 4+i+ £ T . (4) 

As a comparison, the cuts used in the D0 analysis are 

(H t ^t) > (325,225) 2j+ £ T , 

(375, 175) 3j+ # T , (400, 100) 4+i+ ^ T . (5) 

The lowered cuts provide better coverage for intermedi- 
ate mass gluinos, as indicated in Fig. [2] For m g < 200 
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FIG. 2: The 95% gluino-bino exclusion curve for DjZ) at 4 
fb _1 for S/B > 1. The dashed line shows the corresponding 
exclusion region using DjZ) 's non-optimized cuts. The masses 
allowed in the CMSSM are represented by the dotted line; the 
"X" marks the current DjZ) limit on the gluino mass at 2.1 
fb" 1 (see text for details) [7]. The inset shows the effect of 
scaling the production cross section for the case of S/B > 1. 
The solid lines show the exclusion region for a/3 (bottom) 
and 3<7 (top). 



GeV, we place tighter cuts on the monojet and dijet 
samples than D0 does. While D0 technically has 
statistical significance in this region with their existing 
cuts, their signal-to-background ratio is low. Because 
of the admitted difficulties in calculating the Standard 
Model backgrounds, setting exclusions with a low signal- 
to-background should not be done and fortunately can 
be avoided by tightening the Ht and JBt cuts. Similarly, 
for larger gluino masses, the generic cuts are no longer 
effective and it is necessary to use the optimized cuts, 
which arc tighter than D0 's. 



IV. CONCLUSIONS AND OUTLOOK 

In this paper, we describe the sensitivity that D0 
has in searching for gluinos away from the CMSSM hy- 
pothesis in jets + fjT searches. It was assumed that the 
gluino only decayed to two jets and a stable bino. How- 
ever, many variants of this decay are possible and the 
search presented here can be generalized accordingly. 

One might, for example, consider the case where the 
gluino decays dominantly to bottom quarks and heavy 
flavor tagging can be used advantageously. Cascade de- 
cays are another important possibility. Decay chains 
have a significant effect upon the searches because they 
convert missing energy into visible energy. In this case, 
additional parameters, such as the intermediate particle 
masses and the relevant branching ratios, must be con- 
sidered. In the CMSSM, the branching ratio of the gluino 



into the wino is roughly 80%. This is the dominant decay 
affecting the D0 gluino mass bound in CMSSM param- 
eter space (see Fig. 2). While this cascade decay may be 
representative of many models that have gluino-like ob- 
jects, the fixed mass ratio and branching ratio are again 
artifacts of the CMSSM. A more thorough examination 
of cascade decays should be considered. 

In addition to alternate decay routes for the gluino, 
alternate production modes are important when there 
are additional particles that are kinematically accessi- 
ble. In this paper, it was assumed that the squarks are 
kinematically inaccessible at the Tevatron; however, if 
the squarks are accessible, gq and qq production chan- 
nels could lead to additional discovery possibilities. For 
instance, a gluino that is degenerate with the bino could 
be produced with a significantly heavier squark. The 
squark's subsequent cascade decay to the bino will pro- 
duce a great deal of visible energy in the event and may 
be more visible than gluino pair production. 3 

Finally, in the degenerate gluino region, it may be 
beneficial to use a mono-photon search rather than a 
monojet search [22 . 4 Preliminary estimates of the reach 
of the mono-photon search show that it is not as effec- 
tive as the monojet search. This is likely due to the 
absence of final-state photon radiation from the gluinos. 
However, it may be possible to better optimize the mono- 
photon search, because the Standard Model backgrounds 
are easier to understand in this case. 

Ultimately, a model-independent search for jets plus 
missing energy would be ideal. We believe that our ex- 
clusive nj + fjT searches, with results presented in an 
exclusion plot as a function of Ht and fix, would pro- 
vide significant coverage for these alternate channels [21] . 
This analysis should be carried forward to the LHC to 
ensure that the searches discover all possible supersym- 
metric spectra. The general philosophy of parameter- 
izing the kinematics of the decay can be easily carried 
over. The main changes are in redefining the Ht and f^T 
cuts, as well as the hard jet energy scale. We expect a 
similar shape to the sensitivity curve seen in Fig. |] but 
at higher values for the gluino and bino masses. There- 
fore, it is unlikely that there will be a gap in gluino-bino 
masses where neither the Tevatron nor the LHC has sen- 
sitivity. 
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